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Abstract 

The problems of the standard model are briefly reviewed and the motivations for 
introducing supersymmetry are discussed. Two realistic supersymmetric models; the 
Minimal SuperSymmetric Model, MSSM, and its proposed extension NMSSM are in- 
troduced briefly and then some details of the Higgs sector of NMSSM are stated and 
some constraints on parameters are described. I then explore the Higgs boson masses 
and couplings for some interesting scenarios, including a few different ways of taking 
the limit where NMSSM reduces to MSSM. 
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1 Introduction 



The standard model, despite being very successful and in overall very good agreement with 
experiments, has many obvious shortcomings. One of the leading contenders claiming to 
solve some of these problems is supersymmetry (SUSY). As the perhaps most important 
example, we have the problem of the small Higgs mass in the standard model, called the 
hierarchy problem. This problem is automatically solved in supersymmetric theories. From a 
theoretical viewpoint, supersymmetry is attractive not only because it solves some problems, 
but also since it can be said to be the most general possible extension of our normal spacetime 
symmetries [TJ [2] . 

In this study, which is part of a Bachelor thesis at Lund University, some problems 
of the standard model and how they motivate the study of supersymmetry are reviewed. 
Then I briefly introduce the minimal supersymmetric model, MSSM, and one of its proposed 
extensions with a singlet Higgs field, the next-to MSSM or NMSSM. Then I examine the Higgs 
sector of NMSSM in some more detail, and look at some of its qualitative features, using 
numerical tree level calculations to look at the mass spectrum and the reduced couplings. 
Some different limits where NMSSM reduces to MSSM are also studied. 

1.1 Supersymmetry and some problems of the standard model 

The standard model, while highly successful and in good agreement with the absolute ma- 
jority of performed experiments (there are however some persistent measured deviations, like 
the anomalous magnetic moment for the muon[T2"]). still has some serious problems and is 
theoretically very unsatisfactory. For one thing, it has a large number (at least 19) of arbi- 
trary parameters, including the particle masses, the three gauge couplings, the weak mixing 
angle and the CP-violating Kobayashi-Maskawa phase. This is a problem which really isn't 
improved much by supersymmetry (in fact the breaking of supersymmetry seems to introduce 
a lot of new parameters). 

Another obvious deficit of the standard model is that it doesn't describe the fourth known 
force, gravity. What one would ultimately want is a theory that describes all the known forces 
in a unified way. Since supersymmetry is an extension of the ordinary spacetime symmetries 
rather than a new internal symmetry, it is natural how curved spacetime and thus gravity 
follows from making supersymmetry local. The resulting theory is called supergravity, see [3] 
for an introduction. This doesn't really give us a realistic and consistent quantum theory of 
gravity but it at least might hint at how this may be achieved. 

The standard model also has to be modified in order to be consistent with the standard 
theories of cosmology. Especially, it cannot explain the observed cold dark matter, and 
when one calculates the vacuum energy in the standard model it gives contributions that are 
far too large to match the observed small, but nonzero, value of the cosmological constant. 
Supersymmetry offers a good candidate for dark matter, since it turns out that the lightest 
supersymmetric partner (which probably is a mixture of the superpartners of the photon, 
Z-boson and neutral Higgs bosons, called the neutralino) has to be stable, if we assume 
R-parity (which in turn is strongly implied from limits on the proton lifetime). 

Then there is the hierarchy problem concerning the smallness of gravity compared to 
electromagnetism. Since gravity couples to mass, this problem is equivalent to asking why 
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the particle masses are so much smaller than the Planck mass scale (Mp ~ 10 19 GeV), i.e. why 
the ratio of the electroweak scale (~ Mw) and the Planck scale is so tiny, Mw/Mp ~ 1CT 16 . 
From reasons stated above, we know that the standard model is an effective theory, at most 
valid up to the Planck scale. This means that when renormalizing the theory we must have 
a finite energy cutoff, that can at largest be the Planck mass. This finite cutoff is really what 
causes the problem. The masses of the particles come from the Higgs mechanism, so the 
question why the masses of the particles are small is in turn equivalent to asking why the 
Higgs has such a small mass. If one calculates loop corrections to the Higgs propagator, one 
finds that the corrections look like 

bm\ ~ £>(a)A 2 , 

so the correction to the mass squared are proportional to the square of the UV-cutoff A, 
which by our previous reasoning is some large, finite energy, maybe of order Mp or at least of 
some large unification scale (called the GUT-scale, for Grand Unified Theory). So in order to 
give the Higgs its required mass, which of course is much smaller than either of these scales, 
these loop corrections has to be very precisely cancelled. This we can do by giving the tree 
level diagram exactly the required value. Technically, this is not a problem since there is no 
constraint on the value of the bare mass, but it introduces a very heavy fine-tuning into the 
theory, where the bare mass has to have exactly the correct value, and this is theoretically 
very unsatisfying. 

A thing to note is that implicit in this line of reasoning is that you assume that there is 
no need for new physics below some very large energy, i.e. the cutoff scale is large. This is 
called the "big desert" assumption, and while many, holds this to be true, not all physicists 
agree. An argument for this assumption is that in the standard model, the running of the 
couplings is such that at a high energy scale, the grand unified scale (GUT-scale), of about 
~ I0 16 GeV, all the known coupling become roughly the same, which implies that at least at 
this scale, our current physics should drastically change. Supersymmetry actually improves 
this a bit, making the couplings meet more closely than in the standard model. If one 
accepts the big-desert assumption, then the hierarchy problem is real, and some mechanism 
is needed to keep the Higgs mass small. Supersymmetry is the most popular proposal to 
solve this problem, but other theories exists, such as technicolor [U |5] (in which the Higgs is 
a composite particle) and extra dimensions (for example the ADD model[6]). 

A good reason to study supersymmetry is therefore that it solves this hierarchy problem 
in a natural way. It is solved since for every fermion that couples to the Higgs, SUSY 
adds a scalar with the same quantum numbers. When calculating the loop corrections, the 
fermion loops and the scalar loops will be of the opposite sign and (if supersymmetry wasn't 
broken) be of the same size and thus cancel. Even when supersymmetry is broken, these 
cancellations removes the dependence on A 2 and reduces it to a logarithmic divergence [1 lj , 
if supersymmetry is softly broken. 

A final problem of the standard model worth mentioning briefly is the strong CP problem. 
This comes from the observational fact that the strong interaction as described by QCD 
doesn't seem to violate CP symmetry, in contrast with the weak force. This is a problem 
since there are natural terms in the QCD Lagrangian that violates CP conservation. To 
conform to the experimental data, a large amount of fine-tuning is again required. The 
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most well-known proposition for solving this problem introduces new scalar particles, called 
axions[2T], which can make an appearance in the NMSSM, as discussed later. 

Supersymmetry is a vast field, and the theory behind it is an exciting subject. For more 
complete expositions of SUSY and its consequences, there are numerous books and articles 
to consult, for example [TJ, El El EE]- 

2 Realistic supersymmetric models 

In this section I will introduce the simplest way to extend the standard model into a super- 
symmetric theory, the minimal supersymmetric standard model (MSSM). Then an extension 
to this theory is presented, the next-to MSSM (NMSSM), mostly in order to solve a problem 
concerning the value of a dimensionfull parameter in the MSSM. Since nature obviously isn't 
supersymmetric at low energies, we also need to study how supersymmetry is broken. This 
is a large subject which I won't cover in any detail, only introduce the concept of how we can 
introduce so called softly breaking terms into the Lagrangian of our models, as discussed in 
the next section. 

2.1 Softly breaking terms 

The most popular ideas about how supersymmetry breaking works, is that it is spontaneously 
broken, in a manner similar to how the gauge symmetries are broken, see [H] for a discus- 
sion. In fact, supersymmetry is broken as soon as the vacuum gets a nonzero energy. This 
fact means that the breaking of supersymmetry and the breaking of gauge symmetries are 
closely connected subjects. There are different additional terms you can introduce into your 
Lagrangian, such that you can give the fields in these terms a non-vanishing VEV, and thus 
give the vacuum a nonzero energy, breaking supersymmetry. For more details, see |14j . 

When constructing realistic models, we don't really need to care about the details of 
exactly how this happens. Instead, we can introduce terms into our Lagrangian that explicitly 
breaks supersymmetry, but at the same time preserves renormalizability and are such that 
at high energy, above the supersymmetry breaking scale, they become irrelevant. Such terms 
are called softly breaking terms, and are essentially things like scalar mass terms, gaugino 
masses or cubic scalar terms with dimensionfull couplings. For our purposes, analysing the 
Higgs sector at leading order, we only care about the scalar mass terms. So when we have 
the supersymmetric Lagrangian, we can then add all such allowed terms, and view them as 
an effective description of how supersymmetry is broken. 

2.2 The Minimal SuperSymmetric Model 

Just as it sounds, the Minimal SuperSymmetric Model (MSSM) is the model you get when 
you try to minimally extend the standard model to incorporate supersymmetry. Since none 
of the particles in the standard model have the same quantum numbers (excluding mass), 
one cannot let any of the known particles be each others superpartners. So instead we let 
every particle be a part of a corresponding superfield, and then put the superfields in the 
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same SU (2) L doublets as in the SM. The same of course applies to the gauge fields, which 
now become part of gauge superfields. 

In a supersymmetric theory, only a Higgs with hypercharge Y — 1/2 can have the nec- 
essary Yukawa coupling to give masses to the up-type quarks with charge +2/3, and only a 
Higgs with hypercharge —1/2 can have the necessary couplings to give mass to the down-type 
quarks. This is because the superpotential is holomorphic, so the Higgs doublet giving mass 
to the up-type quarks cannot also give mass to the down-type quarks since we are not allowed 
to use the complex conjugate. Thus we at least need two different Higgs S77(2)£-doublets in 
order to give mass to all the massive particles. We will call the Y — 1/2 doublet H u , and 
the Y = —1/2 doublet H d . Then the upper component of the Y — 1/2 will have isospin 
T3 = +1/2 and therefore have electric charge +e. The lower component will be electrically 
neutral, and in the same way the upper component of the Y — — 1/2 doublet will be neutral 
while the lower will have a negative electric charge. So the superdoublets look like 

Ml )> MS) (1) 

and the Higgs fields which gives masses to the fermions will be the corresponding scalar fields. 
We can then note that the Higgs superdoublet H d has the same quantum numbers as left 
handed leptons (and sleptons). Therefore we can use it to give mass to the leptons as well 
as the down-type quarks, so we don't need another Higgs doublet for this purpose. There is 
another way to motivate the need of two different Higgs doublets which is based on anomaly 
cancellations, but this isn't logically needed. 

The minimal superpotential involving these superfields which in a reasonable way extends 
the standard model is 

Wmssm = Vu^iHu ■ Qj - yjdiH d ■ Qj - y^eiH d ■ Lj + fiH u ■ H d (2) 

where the "hatted" letters denote the superfield doublets or singlets corresponding to the 
normal SU(2)l doublets/singlets in the standard model, and i,j are generation indices, 
i = 1,2,3. That is, for the first generation, Qi = (u,d) T , L\ = (e£,z> e ) T , e\ = e.R,u\ = u 
and so on. The Higgs doublets are as described above, and the y 1 - 7 are the Yukawa couplings 
among generations. The products of SU (2) L doublets are given by 

A ■ B = e ab A a B b 

where e a b is the fully antisymmetric symbol in two dimensions with e±2 = 1 and a, b are SU(2)l 
indices. In this superpotential the \x parameter has dimension mass and is contributes to the 
masses of Higgs fields, and it is this simple fact which motivates the introduction of the 
next-to minimal supersymmetric standard model (NMSSM). 

From this superpotential and the ordinary gauge couplings of the standard model, we 
can calculate the scalar potential, by calculating the F and D contributions in the usual way. 
Doing this, and looking only at the Higgs sector of the potential, we find 

V F = fi 2 (\H u \ 2 + \H d \ 2 ) (3) 
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and 

V D = \g 2 W-\H d \^ + \ g l\Hl.H d \\ (4) 

As described above, we can then add the soft supersymmetry breaking terms; 

V soit = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + (m 2 3 H u ■ H d + h.c.) (5) 

where the dimensionfull parameters rn 2 Hu , rn 2 Hd and m\ clearly have to be of the order of 
the weak or supersymmetric breaking scale. The total scalar potential is then the sum of 
these three terms. By letting at least one of m 2 H and m 2 H be negative, H u and H d acquires 
non-zero VEVs, breaking the symmetry. 

From requiring vacuum stability we get some relations between m 2 ^, m 2 ^, m§, the VEVs 
and n; as described in more detail in the next section for the NMSSM. Using these, one can 
calculate and then diagonalize the mass matrices that describe the physical mass eigenstates 
in terms of the parameters of the model. 

In the MSSM, it turns out that we get one physical charged Higgs state, H ± , with a mass 

2 (2m 2 1 2 \ v 2 
m l± = — 1 + -A ~7T (6) 



v u v d 4' 

where (H®) = v u /y/2 and (H d ) = v d /y2, i.e. the VEVs, and v 2 = v 2 +v%, which corresponds 
to the VEV of the Higgs in the standard model. We also get one neutral, pseudoscalar 
(CP-odd) Higgs, called A, with a mass 

where the useful angle (3 is defined from tan/3 = Finally, we also get two neutral scalar 
(CP-even) Higgses, Hi and H 2 (where H\ is lighter than if 2 ), that have the masses 



1 

m H u H 2 - 2 



m A + M z =f \ {m\ + M|) 2 - ^M\m\ cos 2 2/3 



(8) 



If we remember the relation = \{v 2 /2)g 2 Q and express sin 2/3 in terms of v u and v d , we 
see that 

m 2 „± =m 2 A + Mf v . (9) 

We can also conclude that 

m 2 Hl +m 2 H2 = m 2 A + M 2 , (10) 

and the more striking inequality 

m Hl < mm(m A ,M z ), (11) 

meaning that no matter how we choose our parameters, mn 1 < Mz- This is a tree level 
prediction, and loop corrections can lift the mass of Hi above the so far established limits, 
but this is still an important prediction of MSSM. 

1 The extra factor 1/2 coming from the factors of \[2 in my definition of v u and Vd- 
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As we are excluding larger and larger values of m# i; this is really a problem, called the 
little hierarchy problem. Just as the general hierarchy problem, this concerns the separation 
of mass scales, because in order to generate the large loop corrections needed to increase rriHx, 
the other sparticle masses needs to become very large, again creating a new unexplained mass 
scale in the theory. The maximal bound one can get without adding new dynamics to the 
theory is something like mn\ % 135 GeV. Thus results that exclude a Higgs mass lighter 
than that will exclude the whole of MSSM. 

As mentioned above, the /i and mf. parameters are dimensionfull. The mf parameter 
isn't a problem since it enters as one of the softly breaking terms, but \x enters through the 
ordinary Lagrangian, so the only natural values for it before SUSY breaking occurs, is either 
or something similar to the Planck mass Mp. However, to be phenomeno logically viable 
we must have a \x that is of similar size to the electroweak scale. Otherwise there would have 
to be miraculous cancellations between /i 2 and the soft supersymmetry breaking terms. It 
is in order to solve this problem we motivate the study of the NMSSM, where by adding a 
singlet Higgs field and coupling it to the Higgs doublets, the /i parameter is generated through 
supersymmetry breaking. This breaking gives the singlet and the two Higgs doublets VEVs, 
and thus the singlet-doublet-doublet coupling gives us an effective /i. This explains why \x 
should be roughly the same scale as the electroweak breaking scale. 

For a more complete discussion of MSSM, and some discussion about how it may be 
discovered, see for example [TB"] . 



2.3 The Next-to Minimal SuperSymmetric Model 

As stated above, the NMSSM [151 HH] is a proposed extension of the MSSM, which in a 
natural way solves the /x-problem. In order to get rid of the dimensionfull /x parameter, we 
add a new Higgs SU{2) L singlet S to the theory. Of course, in principle nothing forbids a \i 
term just because we add a new singlet, but we take it to have the "natural" value 0. The 
new superpotential looks like 

Wnmssm = Wmssm + ^SH u H d + ^ K ^ 3 , (12) 

where A, k are new, dimensionless parameters of the model. From this superpotential and 
the usual gauge couplings, the F and D part of the potential can be computed. The result 
looks very much like in MSSM, but with some extra terms; 

V F = \XS\ 2 (\H u \ 2 + \H d \ 2 ) + \XH u -H d + KS 2 \ 2 (13) 
Vo = \g\\H u \ 2 -\H d \ 2 Y + \gl\HlH d \ 2 (14) 

Since supersymmetry has to be broken, and we don't know nor care about the details of the 
breaking mechanism, we also have to add to the potential all possible terms which breaks 
supersymmetry in the acceptable, soft way explained in section 2.1 This soft potential looks 
like 

V^soft = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + (\A X SH U ■ H d + ^kA k S 3 + h.c.J . (15) 



7 



Nothing is preventing us from adding the m\H u ■ Hd term present in the MSSM case, but 
this would add an additional parameter with mass dimension, in conflict with the philosophy 
behind NMSSM, so we consider only the case m| = 0. The entire Higgs potential is then 
given by the sum of these, 

Vffiggs = V F + V D + V soft . (16) 

Then, as in the breaking of electroweak symmetry, we assume that m| < so that S = 
is an unstable state. If we define (S) = v s /y/2 we see that when we expand the singlet field 
around its VEV we get a term fi = Xv s /y/2 with mass dimension^} Since this term comes 
from supersymmetry breaking, it's natural for it to have a value of magnitude < MgusY, 
where Mgusy is the scale where supersymmetry is broken. This is the way in which NMSSM 
solves the //-problem of MSSM. 

Further, we also assume that at least one of the other Higgs mass parameters m\ and 
Tn H are negative, so that H u , Hd also get nonzero VEVs, as required to break the electroweak 
symmetry. We then have the gauge freedom to choose = (Hj) = 0, so that the vacuum 
is uncharged. In this treatment, I will discuss the vacuum obtained by further assuming all 
the remaining VEVs to be real, and described by 

(Hu) = 4 ( ° V H d = -±= ( V (S) = ±=v M . (17) 



V2\v u J' a V2 V J ' V2 

We then require this vacuum to be a stable local minimum of the potential, giving us three 
different relations of the type 



dS 



vacuum 



relating the squared masses of the Higgs fields to the VEVs and the other parameters in 
the theory. The derivatives w.r.t. fields with zero VEVs are trivially zero. If solved for the 
masses, these three relations are 

ml = < u + N 2 = ^K-^)+A^(V2^a + ^)-^A 2 (^ + ^) (18) 

m 2 d = m^ + |/i| 2 = ^(^-^) + A|J(v^A A + ^)-^A 2 (^ + ^ (19) 

m 2 s = -\v 2 \ 2 - v 2 s k 2 + -^=A X X^-^ + v u v d n\ 1 —A k v s k. (20) 

^ v2 v s v2 

These new masses (m u , ma) are defined since when we give the singlet a VEV, effectively there 
will be an additional massterm of |/i| 2 for the doublet fields, so we calculate the conditions 
for these effective masses. 

We now see that a full specification of the Higgs sector in the NMSSM requires six 
parameters: A, k, A\, A k , tan/3 and v s . Conventions can be chosen such that A, tan/3 and v s 
(and thus /i) are positive, and this is what I will do. For my purposes I also keep k > 0, since 
switching this sign doesn't change any of my results. In my numerical studies I will replace 
A\ by the physical mass of the charged Higgs, m H ± , which of course must be positive, and 
we will see that the requirement of positive masses squared restrics A K to the negative range. 



2 Note that this definition of /i is a convention, which differs by a factor from the most common one. 
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The "MSSM limit" can be approached smoothly by keeping the ratio k = k/\ fix and 
letting A — > 0, while keeping fx = v s \/\/2 constant. Since the only couplings between the 
Higgs doublets and the new Higgs singlet are dependent on A and k, the singlet field decouples 
in this limit and one recovers the Higgs sector of MSSM. How this works will be explained 
in more detail in section 14.21 

Another thing worth mentioning about the NMSSM, is that it can be used to solve the 
little hierarchy problem of the MSSM. This is because if we let A (or v s ) become larger, 
the mass of the lightest Higgs gets larger, so by having a large A we can get a large Higgs 
mass. This approach is sometimes called A-SUSY[20]. In that model we however give up the 
requirement of perturbativity up to the GUT-scale. 



3 Details of NMSSM 

In this section I will go through some technical details about the Higgs sector of the NMSSM. 
First the mass matrices are described in some detail, and then the couplings of the Higgses to 
the W/Z and the quarks are briefly described, introducing the concept of reduced couplings 
in order to easily compare it with the standard model Higgs and the MSSM. Finally some 
theoretical and experimental limits on the parameter space are discussed. Some other articles 
discussing the Higgs sector of the NMSSM are [TTJ ITS] . 



3.1 The mass matrices 

Since mixing only can occur between states with the same quantum numbers, we get three 
different mass matrices, one for the charged Higgs states, one for the scalar or CP-even 
neutral states, and one for the pseudoscalar or CP-odd states. Since they are obtained by 
taking derivatives of the potential, they are all real and symmetric. This is at tree level, 
taking higher order corrections into account this is no longer the case. The mixing matrices 
being real also means that, at tree level, there is no CP violation. 



3.1.1 Neutral scalar states 

In the natural basis {H u R , H d Rl Sr} where the subscript r denotes the real part of the 
corresponding scalar field, we get the mass-squared matrix as follows 

Ml, = -/vl + ^(V2A x + nv s ) (21) 

M s 2 22 = Ig 2 vl+^(V2A X + Kv s ) (22) 

<ss = v s k^A k + 2v s ^+^\A x V -^ (23) 

Ml 12 = v d v u {\ 2 --g 2 )--v s X(V2A x + v s K) (24) 

M 2 13 = \ (v s (v u \ - v d K) - -^=A x vA (25) 
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M 2 23 = A (v s (v d X-v u K) - -^=A x v u ^j (26) 

where I've used the stability conditions to eliminate the mass parameters from the potential 
in favour of the vacuum expectation values and couplings. This matrix doesn't really lend 
itself to much further algebraic simplification, so it is evaluated in the form given here and 
numerical methods are used to find it's eigenvalues, which corresponds to the physical masses 
of the CP-even Higgs states, which are denoted Hi, H 2 , H 3 , ordered from the lowest mass to 
the highest. 

3.1.2 CP-odd neutral states 

In the natural basis {H® j, H® j, Si} we get the following mass matrix for the pseudo-scalar 
states: 

M 2 n = IV -1 Vs \(V2A x + v s k) (27) 

M 2 22 = h Vs ^y/2A x + v aK ) (28) 

M 2 33 = -1=A k v s k+ V ^\(^=A x + 2kv s ^ (29) 

M 2 12 = ±\v a (V2A x + Kv a ) (30) 
M 2 = lv d X(V2A x - 2v s k) (31) 



2 



<23 



^v u X(V2A x - 2v s k) (32) 



If the first two basis-elements are rotated with the angle (3, a massless Goldstone mode 
decouples, and the new mass matrix (dropping the massless mode) in the basis (P l5 P 2 ) 
becomes 



v 2 

M 2 >n = v s X(V2A x + kv s ) = M\ (33) 

M 2 >22 = V -^(^ + 2v s^)-^A k v s k (34) 

M 2 >12 = v\(A x -2kv s ) (35) 

where we introduce the mass parameter M\. Note that this is not a physical mass, only a 
parameter which can be taken as one of the parameters instead of A x . It can be useful, because 
in the MSSM-limit, Ma becomes the physical mass of the pseudoscalar Higgs. The matrix 
that diagonalises this is of course a 2 x 2 orthogonal matrix, and can thus be parametrized 
by an angle 9a- The new basis in which the mass matrix is diagonal is then 



A x \ _ / cos 9 A sin 9 A \ f Pi 
A 2 ) ~ { -sin9 A cos9 A P 2 



(36) 
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3.1.3 Charged states 

In the natural basis {H^ R , HJ R }, the mass matrix for the charged states looks like 



Ml 12 



vM - 2A 2 ) + ^ (V2A X + kv s ) (37) 
]v d v u {g 2 2 - 2A 2 ) + -Xv s (V2A x + v s k) (38) 



M c 2 22 = \vl(gl/2-2\ 2 ) + V -^(V~2A x + v s K). (39) 

By a rotation through the mixing angle /3, this gives a mass matrix in a new basis {H ± , G ± } 
with the only nonzero element 

v 2 11 

M 2 C , n = v s X(V2A x + v s k) + -g 2 2 v 2 - -X 2 v 2 

2v u v d 4 2 

= M 2 +M^-^AV (40) 

To get the second equality, we use the definition of M\ in addition to the previously noted 
relation (g 2 v/2) 2 = M^. The state G ± is a massless Goldstone mode, which is "eaten" by 
the W ± to give it mass. The charged Higgs state is denoted by H . 

3.2 Reduced couplings 

If we want to express how the physical Higgs particles, i.e. the mass eigenstates, couple to 
fermions and gauge bosons, what one needs to do is to express the original weak eigenstates 
H u , Ha, S in terms of the mass eigenstates H ± , A%, A%, Hi, H 2 , H 3 . This is of course done by 
looking at the matrices that rotates the weak eigenstates into the mass eigenstates, i.e. the 
mixing matrices as defined above. 

We are primarily interested in how the \^-boson couples to the different Higgses, where 

V can be either W ± or Z, and the couplings to quarks, since in the generic Higgs decay 

2 

H //, there is a factor -^f- meaning that the heaviest fermion allowed dominates, i.e. 

m w 

either the top or bottom quark. 

The way to find these couplings is to write down the relevant terms in the Lagrangian, 
which is originally in terms of the weak eigenstates, and then re-express it in terms of the 
mass eigenstates H, H ± and Ay The details can be found in [T8] . 

In order to simplify the notation and keep it from getting unneedingly cluttered, we define 
so-called reduced couplings, where we take the full coupling and divide out the associated SM 
coupling, 

n _ dVVHi n _ gZAiHj _ gZA z H 3 

^vvh, = -snT> <^za z h 3 - rSM - , , 

UVVH ^ZHH y/ z 

where g = \Jg\ + g\, gi and g 2 being the gauge couplings of the electroweak force. Here, 

V can stand for either W or Z, the reduced coupling will be the same in either case. The 
reduced coupling we will look the most at is the H^VV coupling, since this measueres how 
standard model like the scalar Higgses are. If we let Si = (H uR , H d R , Sr) be the scalar weak 
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eigenstates, and Hi = J2j (i- e - &ij is the mixing matrix), then this reduced coupling is 

defined as[TT] 

G^vv = sin PSh + cos (3S i2 . (41) 
The H i H ± W T coupling is similarly given by 

G Hi h±wt = cos /3S a - sin /3S i2 . (42) 

Since these reduced couplings comes directly from orthogonal mixing matrices, we may 
conclude that they should fulfil certain sum rules. This is because of the sum rules that 
elements of orthogonal matrices fulfil: the sum of the squares of one row (or column) is equal 
to one. For the reduced couplings, this means 

i i 

In the same way, the reduced couplings of the Higgses to the top and bottom quarks also 
comes directly from an orthogonal matrix, but in this case a dependence on tan /3 also enters, 
since this describes how large the difference is between the two VEVs of the Higgs doublets. 
In this case the sum rules are 

y] G 2 H tt = . 2 , G 2 Hlbb = 2 „ ■ 

z — ' sin a L — 4 cos z p 

i i 

We also have the sumrules from the columns, for example 

G 2 HiW + G 2 HiH±w^ + ^3 = 1j i = 1, 2, 3 (45) 

where S i3 is the singlet component of Hi. If S i3 ~ 0, then Hi will be purely doublet and 
the corresponding sum rule G 2 HiVV + G 2 h . h ± W t = 1 is recovered. Conversely, if Si% ~ 1 
then both the other couplings will be suppressed, which means that detection of Hi will 
be difficult. These sum rules are quite trivial in nature, but can be a useful check on the 
numerical methods used. They are also important phenomenologically, since they in effect is 
a good measure how standard model like the different Hi,Ai are. For example, in this last 
sum rule, if the first term is large it means that Hi is SM like, if the second term is large 
it means there's a large coupling between the doublets making Hi MSSM like, and the last 
term corresponds to how singlet-like Hi is. 

Of course, since the standard model only has one (scalar) Higgs particle, only the H 1 
couplings (assuming that the lightest Higgs also will be the standard model like) have a 
direct correspondence in the standard model. Nevertheless we can define reduced couplings 
by scaling away the gauge couplings and masses. 

3.3 Constraints on the parameters 

In this section I will explain some theoretical and experimental limits on the parameter space, 
and motivate the choices of parameters later used when studying some numerical results. I 
will discuss for which intervals it is sensible to choose values for the parameters, which I 
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choose as A, m#±, k, tan and A K . Some limits can be found from theoretical considerations 
and requirements, while others come from experiments at accelerators or astrophysics. 

As said above, I use the value v = 246 GeV for the electroweak scale, and choose the 
'standard' value /x = Xv s /y/2 = 200 GeV, which allows me to see v s as a fixed value when 
A is chosen. There is a restriction from LEP [T9] on the minimal size in /i, requiring that 

> 100 GeV, coming from lower limits on Higgsino masses, but this is really for MSSM. 
Nevertheless, a too small value of doesn't work. 

Since my analysis is at tree level, we will not discuss parameters entering at loop level, 
where the Higgs masses get corrections depending on for example the top and the stop masses. 
There is probably many cases where even at tree level limits from measurements could be 
used to rule out large parts of the parameter space, but doing this in detail is regretfully 
beyond the limited scope of this study. 

A general way to restrict the parameter space is to require that they fulfil some grand 
unified scenario where all couplings of the same type gets the same value at the GUT scale. 
This is called universal boundary conditions, and will not be required here. 



3.3.1 A and k 



We can see that with k = 0, the Lagrangian, (12), has an additional U(l) symmetry, called 



Peccei-Quinn symmetry [2T] (henceforth called PQ-symmetry). This symmetry was proposed 
as a solution to the strong CP problem, i.e. the problem of explaining why QCD doesn't 
seem to violate CP symmetry like the electroweak interactions do. If this symmetry is exact, 
i.e. k — 0, it will be spontaneously broken by the nonzero VEV of the singlet scalar, which 
will give rise to a massless Goldstone boson, called the Peccei-Quinn axion. This axion will 
show up as the extra pseudoscalar Higgs field (compared to in the MSSM). However, this 
case can in principle be ruled out since it would since there are lower bounds on allowed 
axionmasses[12j which only can be avoided if 10 -16 < A < 10~ 7 [18j. Such a small value of 
A would mean that v s would have to grow very large, making the model unattractive as a 
solution to the //-problem. 

So from this we conclude that we need a nonzero value of k, breaking the PQ symmetry. 
The size of the k coupling will regulate how badly this symmetry is broken, and with a 
small value, only slightly breaking PQ symmetry, we will get a nonzero mass for the lightest 
pseudoscalar. 

If one uses the requirement that A, k and the Yukawa couplings should stay small (eg. < 1) 
so that perturbation theory can be used up to the GUT scale, and uses the renormalization 
group flow, one can get the approximate limit at the electroweak scale |18j 



VA 2 + k 2 < 0.7 (46) 

Also, from choosing a large number of different values of A and k at the GUT scale and using 
the renormalization group equations to run them down to the electroweak scale, one can see 
that the flow favours a small k value. 

We also note that if A gets too small, this forces v s to become big, which means that the 
model no longer works well as a solution to the //-problem. Even if we allow v s to take a 
value of a few TeV, say 2 TeV, which is well over but still "close" to the electroweak scale in 
some sense, this places a limit on A > 0.1, so we get a rather stringent condition on A. 



13 



> 

CD 

o 



05 
CO 




-600 -400 -200 
A_kappa [GeV] 



200 



Figure 1: The masses as a function of A K , where tan/3 = 7, m#± = 250 GeV, k 
X = 0.3. 



0.3 and 



If universal boundary conditions at the GUT scale are imposed (which gives us the so 
called constrained NMSSM|22J), we also get that the ratio X/k has to be close to 3. This is 
however not something that will be exclusively used since I don't in general impose universal 
boundary conditions. 

3.3.2 tan/3 and A K 

The range of A K is rather tightly constrained from the condition of vacuum stability. In figure 
[TJ the only allowed range of A K is where all the masses are positive, i.e. —600 < A K < —30 
GeV. In many cases the limits are a lot stricter than this. From such plots you can also 
see that for some choices of the other parameters there are no acceptable value of A K at 
all; in some cases the lightest pseudoscalar and the lightest scalar never both get a positive 
mass at the same time. For the coming plots where the A K dependence matters, the value 
of A K is chosen roughly in the middle of its allowed range, for a typical value of the running 
parameter. The typical values are —100 and —250 GeV. 

As for tan /3, an analysis of the running couplings shows that a low value of tan /3 is 
favoured. However, experiments rule out a too small value, so a not so small value is required 

P2J. 

The theoretical upper bound is tan/3 < m t /m b ~ 50, and I will briefly study what 
happens when you take a large value, tan/3 = 30 in the model. See also [23] for a study of 
what happens when you saturate this upper bound. We also have a lower bound tan/3 > 1.2 
from requiring At, the top quark Yukawa coupling, to remain small up to the GUT-scale [1 7] . 

3.3.3 m H ± and the other Higgs masses 

For the Higgs masses experiments have placed general lower limits. For the neutral Higgses, 
LEP has published negative search results [19] in some different decay channels, and depending 
on the precise branching ratios the limits looks a little bit different, but generally the lower 
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bound from LEP is around m Ht > 80 - 90 GeV, for the MSSM. In the NMSSM these limits 
can be avoided, but they still give some kind of general guidelines. In the numerical studies, 
a lower value of = 90 GeV, and a higher value of 250 GeV will be used, when we don't 
let m H ± vary. 

From different experiments at the Tevatron and LEP, rather strict limits can be placed 
on a SM-like Higgs [T2]. for example the Tevatron has excluded the range 158 GeV < rrih < 
173 GeV, and from LEP we have exclusions of a mass lower than 114 GeV, but since the 
couplings to fermions and gauge bosons of the MSSM or NMSSM Higgs bosons can be 
suppressed compared to the SM Higgs, these limits can be avoided. 

4 Results 

In this section I will present numerical results that explores some of the features of NMSSMs 
parameter space. First, masses and couplings are treated as a function of the charged Higgs 
mass. From this we see some possibly interesting features of the model. Then it is studied 
how the masses and reduced couplings (and thus the mixing) behave in some different kinds 
of MSSM limits. 

Since my calculation is only at tree level, it is not sensible to compare directly with 
experimental limits. Even so, the general features are maybe even better understood at tree 
level, since it is easier to compare directly with the formulae without too much cluttering of 
the expressions. For my numerical results, I've written code in Java, using the basic linear 
algebra library Jama to diagonalize and find eigenvalues of matrices. 

4.1 Varying the charged Higgs mass 

In order to see how the NMSSM mass spectra behaves, it can be instructive to plot the masses 
as a function of the charged Higgs mass. From these plots, and the requirement of vacuum 
stability (i.e. m 2 H > 0, the lightest scalar mass positive) we can find limits on allowed values 
for mn± for fixed values of the other parameters. 

As an aside, just in order to confirm the theory and my numeric calculation, we can check 
that the sum rules for the reduced couplings actually holds in practise, which it turns out 
they do. A thing to note when looking at plots of reduced couplings and cos 9 a is that what 
I really plot is the absolute value of the couplings. This is for two reasons, first of all that all 
we really care about is the strength of the different couplings, the sign can of course matter 
(mostly when you go to higher orders) but not for our purposes here. The other reason is that 
the numerical method used switched signs discontinuously, so without taking the absolute 
value the graphs looks very discontinuous and strange. This can probably be fixed rather 
easily, but since the sign doesn't matter for our purpose no effort was expended on this. 

4.1.1 The NMSSM with a small « 

In figure[2]we can see that for the following choices of parameters, A = 0.3, k — 0.1, tan (3 = 2, 
A K = —100 GeV, m H ± has to be between 360 and 550 GeV. This case is representative for 
a small K-value, which corresponds to a slightly broken PQ-symmetry. This is favoured by 
the renormalisation group flow. From the figure we also see that the mass of the lightest 
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Figure 2: The Higgs masses, some different scalar couplings and cos 8 a as functions of m^±, 
for k = 0.1, A = 0.3, A K = —100 and tan/3 = 2. In (b,c,d), the physical range of the parameter 
space is inside the two black lines. is the singlet component of the H t state. 
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pseudoscalar almost doesn't change as soon as m H ± gets above ~ 250 GeV. This is the singlet 
dominated pseudoscalar whose mass mainly comes from the ^kS 3 term in the Lagrangian. 
The heavier states that grows more or less linearly are the fields with little mixing with the 
singlet, (as seen from that S33 and cos #,4 both are small) and behave the same way as in the 
MSSM. For these parameters, the next lightest scalar state, m H2 also doesn't depend very 
strongly on m#±, and grows much slower than in MSSM. 

So we see that even though the spectra of high mass states stay roughly the same, the 
three light Higgs states means that a NMSSM with parameters close to these will be easily 
distinguished from MSSM, even if we only find the lighter Higgses. Of course, this only works 
if the reduced couplings, Gh.vv an d Ga^z doesn't become too small to prevent detection, 
which happens in a small part of the relevant parameter space. 

In figure [2}b the couplings of the scalar fields to the W/Z bosons as well as cos6U are 
plotted as a function of the charged Higgs mass. Since the reduced couplings can be thought 
of as a measure of the mixing between weak eigenstates, this plot shows that the mixing 
depends on the charged Higgs mass in a slightly complicated way. The lightest Higgs is as 
one might guess the most standard model like, and it is the mixing with H 2 which raises 
m 2 H above zero. For a specific value, m H ± = 466 GeV, the coupling Gh 2 vv = 0, so if this 
specific scenario is true, the next-lightest Higgs would be totally singlet-like (as we see in 
figure [2jd). In this case, the H 2 Higgs would be totally undetectable through the channels 
used to look for the standard model Higgs. We also see that in the physical range at least, 
the pseudoscalar mixing only varies a little. The H i H ± W T couplings, which measure how 
doublet or MSSM-like the scalars are, vary very little in the physical range, but we do see 
that H\H ± W T pass through zero when m#± = 449 GeV. We also see that the heavy H% is 
doublet-dominated. 

From the figure we can also see that tua 1 and seems to switch behaviour with respect 
to the charged Higgs mass around m H ± ~ 150 GeV. This switch is also apparent in how cos 9 a 
behaves. After this however, the pseudoscalar mixing stays more or less constant, and doesn't 
vary rapidly in the physical region as the scalar mixing does. In the same way we see the 
switch in behaviour between m# 2 and rriH 3 reflected in how all of the corresponding couplings 
switch, although this happens outside the physical region. 

4.1.2 Larger k 

If we let the value of k get larger, the PQ symmetry is more badly broken and the lighter 
pseudoscalar gets a larger mass. This is not favoured by the renormalization group flow, but 
we have no a priori reason to exclude it. In figure [3] we have plotted the mass spectrum 
and couplings as functions of for k = 0.5. In this case, the lightest Higgs is the most 
standard model like by far, and H 2 , H 3 again switch behaviour, around m#± ~ 520 GeV. 

We also see that this large k loosens the constraints on m H ± from vacuum stability. The 
value of A K = —500 GeV used in the figure has been chosen approximately in the middle of 
its allowed range for these parameters. 

In this case, compared to the previous case with k small, we see that apart from the lightest 
Higgs, the rest of the masses are significantly larger. However, they are not extremely heavy 
and are still very much within the range of detection, but the spectrum of light Higgses present 
in the previous case is absent. This would make it harder in this case to distinguish between 
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Figure 3: The Higgs masses, some of the scalar couplings and cos 9 a, as functions of m H ±, 
for k = 0.5, A = 0.3, A K = —500 and tan (3 = 2. In (b,c,d) the physical range is again between 
the black lines. 
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Figure 4: The Higgs masses and reduced couplings of the scalar states to the Z and W bosons, 
as functions of m#±, for a higher value of v s , equivalent to \x = 1000 GeV, instead of the 
normally adopted value of /i = 200 GeV. The rest of the parameters are the same as in figure 
[2j In (b), the physical range is inside the vertical black lines. 



the MSSM and the NMSSM if only the two lightest Higgses can be detected, compared to 
the case with smaller k. 

Just as in the previous case, the pseudoscalar mixing, i.e. cos^a, stays rather constant 
except at the place where the two pseudoscalar fields switch identity. In this case this 
switching behaviour is more distinct, something we see both in how the masses and how 
cos 9 a behaves. 



4.1.3 A larger v s or /i- value 

If we vary the expectation value of the singlet field, v s , or equivalently the value of the effective 
H = v s X parameter, this doesn't change the qualitative behaviour of the mass spectrum very 
much, but it changes the quantitative behaviour. All but the lightest Higgs gets heavier, 
including the charged Higgs, since the region of vacuum stability gets pushed upwards, see 
figure 4 We also see that the constraint from vacuum stability is relaxed (see for comparison 
figure 2), and that the charged Higgs mass, as well as the masses of if 3 and A2 are now forced 
to be larger than ~ 1.5 TeV. Since the lightest of the Higgses remains light, it means that a 
higher v s must make the Hi more SM-like, which means that the H2 becomes more singlet- 
like. That Hi becomes doublet-like means that its coupling to W and Z should become large, 
and this is indeed also the case, as we can see in the right panel of figure |4j 

If we make the same plots for a large value of k, the effect on the mass spectrum will be 
bigger and the heavy states will thus be even heavier, since it gets amplified by the large k 
value. In this case the heavy singlet dominated fields will decouple and the lower mass states 
will behave like in the MSSM, so in this case the distinction between NMSSM and MSSM 
will be hard to find. 
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Figure 5: The Higgs masses as a function of for tan f3 = 35. In (a), we have k 
0.1, while in (b) k = 0.5, A = 0.5. In both cases A K = -100 GeV. 
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4.1.4 Large tan (3 values 

What happens if you increase the value of tan /3, making one of the doublet VEVs much 
larger than the other? In figure |5j this is shown as a function of the charged Higgs mass for 
two different choices of the singlet couplings. From the figure we see that there seems to be 
rather distinct points where the identity of two different Higgses seem to switch. For example, 
the Hi and the H 2 states seem to switch behaviour with respect to the charged Higgs mass 
at around m H ± = 120 GeV. This kind of behaviour really comes from the way the mixing 
matrices depends on m#± (or equivalently A\) and from how we label the different states. A 
large tan (5 value means that the mixing with the H u doublet will be much more important 
in terms of mass than the mixing with the Hd doublet. From this we can understand why a 
larger tan /3 value makes the identities of the Higgses more sharply defined. 

What is perhaps more interesting to note is that we have three masses here that are almost 
independent of the charged Higgs mass, even though exactly what we call the state varies 
with m,H±- This is also coupled to the fact that the couplings in this plot is rather small, 
k = A = 0.1. If we instead make them larger, we instead get the behaviour seen in figure |5j 
where the "switching" behaviour is not at all as sharp. Increasing the singlet couplings also 
pushes the lowest physically allowed value for m H ± upwards (i.e. the value where m Hl > 0), 
and seems to give the lightest scalar Higgs a very small mass. 



4.2 The MSSM limit 



A few things can be noted analytically when we take the limit A, k — > while keeping k/\ = k 

4 



fixed. For example, we can see that the parameter m| in the Lagrangian will approach a 



fixed value as soon as A and k get small. This is seen by looking at the expression for m| we 
got from the requirement of vacuum stability, equation (20). Since we keep \i = ^v s \ fixed, 
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Figure 6: How the masses vary when A goes to zero while keeping k fixed. In this plot, 



tan (3 = 2, m H ± = 250 GeV, k = 0.1 and A h 



-100 GeV. 



we have that 



rrit 



2 \ 2 22, 1 a \ V u v d , 1 A 
-V A —V„K H -=A\A \-V u VdKA =A K V S K 

2 s s/2 V s y/2 



— > kfi(A K — kjj) as A — > 0, 



also using the relation between k and A. If we instead keep k fixed and let A — > 0, we see 
that instead 



1 2 \ 2 2 2 i a \ 2 v uVd \ A ^ 
-V A -V S K + -A X X h V u VdKA — A k k — 

2 2 fi A 



—7- oo, as A — > 0. 

This behaviour is shown in figure [6j and it means that just as we can see from looking 



at the Lagrangian (12), as A goes to zero the singlet field decouples and the masses of the 
singlet dominated states blow up. 

Another way of taking an interesting limit is to let k — > keeping A constant, as can be 
seen in figure [7j Here we see that in this limit (which really isn't a proper MSSM limit since 
A stays large and thus the singlet doesn't fully decouple), the Ai state becomes massless, 
which again is because we restore the PQ symmetry turning A\ into the massless PQ-axion. 
However, it is seen that for all cases with fixed A there is no way of keeping H\ at a positive 
mass squared as k goes to zero. For smaller values of A, m 2 Hi becomes negative for smaller 
values of n, but any given value of A ultimately restricts the lowest possible value of k. So 
from looking at this in addition to the above discussed A — > limit we see that in order 
to decouple the singlet and reduce the theory to MSSM, one is in effect forced to take the 
simultaneous limit k, A — > (or in addition take the limit A R — > as studied later). 
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Figure 8: The masses, i^V^-couplings and cos6U as the MSSM limit is approached, using 
k = 1, tan/3 = 2, m H ± = 250 GeV and A K = -100. 



22 




0.01 

lambda 

(a) Masses 





1 .2 ■ 




1 ■ 


CO 




□> 






0.8 ■ 










o 




o 

X! 


0.6 ■ 


0} 




O 




13 
X) 


0.4 ■ 


en 


CC 






0.2 ; 




- 



H1VV 
H2VV 
H3VV 
cosThetaA 




0.001 0.01 0.1 

lambda 

(b) Couplings and cos 8a 



Figure 9: How the masses, cos 9 a and HiVV couplings vary when approaching the MSSM 
limit, using k = 0.2, tan p = 2,m H ± = 250 GeV and A K = -100 GeV. 



In figure |8j we see that when approaching the MSSM limit in the sensible way, none of 
the masses becomes large, which is understandable since we above showed that m| tends 
towards a finite value. However, also in this case the singlet decouples and does not mix 
with the other fields, which we see by looking at the reduced couplings and cos^a, seeing 
that the lightest Higgs state becomes completely standard model like, and also cos6U goes 
to 1, so that both mixing matrices become block diagonal and there is no mixing between 
the doublet and the singlet. That the lightest scalar becomes standard model like in this 
limit means that for small k and A of roughly the same size (k ~ 1) means that it could be 
detected as easily as in the standard model. 

We also see that cos 9 a goes to 1 much slower than the reduced scalar couplings. This is 
not a general feature but depends on the value of k. However, I've not found any cases with 
k > 1 where the pseudoscalar mixing disappears slower than the scalar mixing. This means 
that small values of A and k suppresses mixing between scalar singlet and doublet states 
much more than between the pseudoscalar states. Moreover, as A becomes small the mixing 
between the scalar Higgs fields stop depending on A and stays more or less constant. This is 
of course very reasonable since A determines the coupling strength between the singlet and 
the doublets. 

However, as we can see in figure [9j the reduced couplings (and thus the mixing) doesn't 
always vary slowly when we take the limit. This behaviour seems to occur only when k < 0.5, 
in the figure we have k = 0.2 as a representative case. In these scenarios, the couplings 
continue to vary very rapidly (considering the logarithmic scale) even when A is very small, 
and the mixing only disappears when A becomes really small. Differently from before, the 
lightest Higgs isn't the most standard-model like in this scenario. This role is instead filled 
by H 2 . We also note that the masses of Hi and H 2 gets very similar in the MSSM limit. 

Since cos 9 a approaches 0, we see that the lightest pseudoscalar state also decouples 
in the MSSM limit (that cos6U = of course means that the off-diagonal element of the 
pseudoscalar mixing matrix sin 9 a = 1). 
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Figure 10: How the masses vary when approaching the MSSM limit for two different values 
of k. The rest of the parameters have the same value as those in figure |8| 



We however see that the masses of the two lightest Higgses are very small for these 
parameter choices, so this particular case is not realistic. For such small values of k the 
requirement that m 2 Hi > seems to rule out most of the parameter space, i.e. for many 
choices of other parameters, m 2 H < 0. I've not found any case where the same thing happens 
for the pseudoscalar mixing, but no methodical search of such a scenario was carried out. 

As we can see from the formula for m| in the MSSM limit, if we let the ratio k = k/\ 
get larger, then the masses of the singlet states should increase. And this is exactly what 
happens as we can see in figure 10, where a larger k value is seen to push the heaviest scalar, 
the singlet dominated one, far up in mass, whereas the mass of the heavy pseudoscalar is 
also increased but not at all as much. 

Another possible way of approaching the MSSM limit is to also send the A K parameter to 
zero, which effectively further suppress the S 3 term in the Lagrangian and thus restores the 
PQ-symmetry. And as we can see in figure 11 in this case the lightest pseudoscalar Higgs 
become massless in the limit, restoring the massless PQ axion. In this limit, we also see that 
the pseudoscalar mixing disappears (cos 6 a — > 1) and that the lightest scalar Higgs again is 
fully standard model like. 

One could also imagine letting v s go to zero smoothly, but in order to keep /i in an 
acceptable range this would mean that A would have to become larger than allowed by the 
requirement of perturbation theory being valid up to the GUT-scale, i.e. A < 0.7 as stated 
before. A small value of \x is excluded from bounds on the chargino masses [12] from direct 



searches. As said towards the end of section 2.3 , we can give up the perturbative requirement 
and consider A ^> 1 as in the A-SUSY model [20], in which case this limit could be viable. 



5 Summary and conclusions 

In this study, we have briefly reviewed some of the motivation behind studying supersymme- 
try. Then the MSSM was described, and its simplest extension, the NMSSM was introduced 
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Figure 11: How the masses, couplings and cos#a vary when A goes to zero while keeping 
k = n/X fixed, in addition to also sending A K to zero. In this plot, k = l,tan/3 = 2,m#± = 
250 GeV. A K begins on a value of —250 GeV, and is sent to zero in the same manner as A 
and k (i.e. we keep the ratio n/A K constant). 



as a solution to the /i-problem. Some technical details of the Higgs sector of NMSSM was 
stated, including the mass matrices and definition of the reduced couplings. This was followed 
by some numerical studies of different choices of parameters, including a look at approaching 
the MSSM limit in some different ways. 

We find that in the by renormalization group flow favoured choice of k, tan/3 and A, the 
mass spectra with three different relatively light Higgs bosons should make the theory easy 
to distinguish from the MSSM even when not all the Higgses are detected. But in other 
perfectly allowed cases, the distinction might not be directly obvious. We also see that for 
this case some of the couplings to SM-particles pass through 0, so that it is possible that for 
example the H 2 state can be hidden and not interact in a standard model like way at all. 

For a larger k value, i.e. a more strongly broken PQ-symmetry, we see that the couplings 
behave in a qualitatively different way. In this case the switching behaviour takes place inside 
the physically sensible area, but the couplings of the Hi state doesn't show any complicated 
dependence on m#±. Over the whole range, Hi is the essentially fully standard model like. 

If we want to study the limits where the singlet decouple, from looking at what happens 
when only one of k or A are sent to zero, we conclude that sensible limits exists only when 
both of them are decreased simultaneously. When we approach this MSSM limit, keeping 
the ratio k/X constant, the mixing with the singlet field disappears. Mostly this decoupling 
happens rather quickly, but in some cases, when k/X < 0.5, the mixing of the scalar states 
depends strongly on the nonzero A. We also see that in this case, the lightest Higgs is no 
longer necessarily the most standard model like, that role being taken by H 2 . But for A > 0, 
both Hi and H 3 are slightly standard model like, so if this is the case (i.e. we have a small A 
and a ratio as described), we could detect three relatively light Higgs bosons with different 
masses. However, these cases seem to depend very much on the ratio k/X having a specific 
value, and also seem to give the lightest states too low masses for it to be realistic, but 
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since I have not scanned all of parameter space and in addition I am only doing tree level 
calculations, this kind of scenario cannot be altogether ruled out. 

So we see from all this that the NMSSM model offers many interesting possibilities not 
seen in the MSSM. 
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